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ABSTRACT 

We present final results on the angular power spectrum of total intensity 
anisotropies in the Microwave Background from the Cosmic Background Imager 
(CBI). Our analysis includes all primordial anisotropy data collected between 
January 2000 and April 2005, and benefits significantly from an improved max- 
imum likelihood analysis pipeline. It also includes results from a 30 GHz fore- 
ground survey conducted with the Green Bank Telescope (GBT) which places 
significant constraints on the possible contamination due to foreground point 
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sources. We improve on previous CBI results by about a factor of two in the 
damping tail. These data confirm, at ~ 3a, the existence of an excess of power 
over intrinsic CMB anisotropy on small angular scales {£ > 1800). Using the 
GBT survey, we find currently known radio source populations are not capable 
of generating the power; a new population of faint sources with steeply rising 
spectral indices would be required to explain the excess with sources. Exten- 
sive testing does not reveal any instrumental effect capable of giving rise to 
the observed excess. We also present a full cosmological parameter analysis of 
the new CBI power spectrum, the WMAP 5-year data, and the latest ACBAR 
data, self-consistently including in the analysis the foreground signal from the 
Sunyaev-Zel'dovich Effect (SZE) from galaxy clusters. With CBI alone, the full 
parameter analysis finds the excess is 1.6a above the level expected for a (Xg = 0.8 
universe. We fit two different SZ templates to the power spectrum and find they 
give markedly different inferred as values suggesting more theoretical work is 
required. We find the addition of high-£ CMB data substantially improves con- 
straints on cosmic string contributions to the TT power spectrum as well as the 
running of the scalar spectral index rirun, but that rirun is quite sensitive to the 
level of the SZE. We also present forecasts for what other experiments should see 
at different frequencies and angular resolutions given the excess power observed 
by CBI. We find that the reported high i bandpowers from current high resolu- 
tion CMB bolometer experiments are consistent with each other and CBI if the 
excess power is due to the SZE at the CBI-level of 2.5 ± 1 times the ag = 0.8 
standard SZ template. This is not the case if the CBI excess source has a flat 
frequency dependence in thermodynamic temperature. 

Subject headings: cosmology, cosmic microwave background 



Introduction 



Measurements of Cosmic Microwave Background (CMB) anisotropies in total inten- 
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confirmation of tlie inflationary structure formation paradigm. Key elements of tliis pic- 
ture wfiicli liave been confirmed are: tliat tlie universe is spatially flat; that anisotropies in 
the microwave background formed via simple causal processes from an approximately scale- 
invariant spectrum of (probably adiabatic) primordial inhomogeneities; and that the grav- 
itational instability picture of subsequent structure formation from the collapse of baryons 
and dark matter in an expanding metric. In spite of what must, on the whole, be de- 
scribed as stunning experimental confirmation of the inflationary predictions, several sur- 
prises emerged, and tensions and ambiguities in our understanding of the data remain. Chief 
amongst these are the necessity for a (presently dynamically dominant) dark energy compo- 
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2008). A number 



anisotropies — pre- 
dominantly the Sunyaev-Zel'dovich Effect (SZE) from galaxy clusters — are expected to 
contribute on small angular scales, thus providing a view of the more recent evolution of 
large-scale structure through the CMB. 

The Cosmic Background Imager (CBI) is a 30 GHz interferometer that measures CMB 
an isotropies from ^ ~ 400 to £ ~ 3000. The instrument has been described in detail 

From January 2000 through No vember of 2001 the CBI surveyed 
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ily on CMB polarization observations, which were conducted from September 2002 through 
April 2005. The fields observed in the polarization campaign encompassed ~ llSdeg^ in all, 
and partially overlapped with the fields observed in the total-intensity campaign. The total 
area covered in the combined datasets is 143deg^. While roughly half of the CBI baselines 
were cross-polarized for these observations (appropriate to measure polarization), the other 



half were co-polar and thus improve r esults on the TT 



polari zation campaign are presented in lReadhead et al.l (l2004bl ) (Paper 8) and lSievers et al. 



power 



spectrum. Res ults from the 



(120071 ) (Paper 9). In this paper we combine all of these data to present the final TT power 
spectrum from five years of CBI measurements. We also examine the implications of the 
observed high-£ signal in the context of a full cosmological parameter analysis including a 
treatment of uncertainties in the SZ foreground models. 

A key limiting factor in interpreting the high-£ excess measured by CBI has been un- 
certainties associated with the extragalactic point source correction. Extragalactic sources 
reduced the sensitivity of the CBI at high-£ in two ways: by requiring a substantial amount of 
data to be "thrown out" owing to possible contamination from known sources, and through 
the uncertainty in the power spectrum of fainter sources extrapolated from number counts. 
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In order to address the second issue, we conducted a 30 GHz survey with th e GET and the 
OVRO 40-m telescope, covering a total of 3562 NVSS (ICondon et al.l 119981) sources in the 
CBI fields; this survey is described in the companion paper iMason et al.l (120081 ) . This survey 
has allowed us to place much tighter constraints on the point source foreground contribution 
to our power spectrum, which we have quantified and included in this analysis. Because of 
the potential of source variabihty, we do not attempt to reclaim the sky area under NVSS 
sources that are not detected by the GET. 

The structure of the paper is as follows. In § [2] we briefiy summarize the data. § [3] 
provides a description of improvements to our analysis algorithms, describes tests performed 
on the data, and summarizes our knowledge of foregrounds (chiefiy discrete sources) in 
the CEI fields. In § H] we present the CEI power spectrum, discuss the significance and 
characteristics of the high-£ excess signal, and give constraints on cosmological parameters 
from this. § [5] summarizes our findings. We describe our spectrum-fitting procedure in 
Appendix |Al 



2. Observations 

For this analysis w e use data collected on both the CEI total intensity fields from 



m 



Readhead et al.l (l2004af) and the intensity data from the CEI polarization fields described 



Sievers et al.l (120071 ). The fields, observing strategies, data reduction, and calibration are 
described in detail in Papers 2, 3, 7, 8, and 9; here we simply summarize the contents of 
these papers. 

The CEI total intensity and polarization fields were spaced by ~ 6 hours in Right As- 
cension (centered near 02h, 08h, 14h, and 20h) and are near the celestial equator. The 
polarization and total intensity fields overlap but are not identical. To test for systematics 
and obtain higher signal-to-noise ratio at high-£, some very deep integrations on individual 
pointings were performed in each campaign. During the total intensity campaign, all in- 
tegration time at 8h was concentrated in a single pointing; the 14h and 20h mosaic fields 
also contained one pointing each on which significantly more integration time was acquired. 
During the polarization campaign, observations at 20h concentrated on a single 1x6 strip 
instead of 6x6 mosaics. These fields are summarized in Table [H The deep 14h and 20h 
fields, while within the total intensity mosaics of these regions and collected and analyzed 
identically, are listed separately in this table since for some tests we analyze only the deep 
data. 



The largest source of systematic errors in the CEI data is a ground spillover signal, which 
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is strongest on the short basehnes. We find this signal is stable over timescales of at least 20 
minutes. For the total intensity observations, contaminating ground signal was removed by 
differencing pairs of individual pointings separated by 8 minutes in Right Ascension. For the 
polarization mosaics 6 successive pointings were observed in succession and a single common 
mode removed from matched visibilities. 



All data are calibrated with respect to the 5- year WMAP Ka Jupiter brightness temper 



atur^ Tjupiter = 146.6 ± 0.75 K (IHill et al.ll2008l ). We adopt a 1% in amplitude calibration 

147.3 ± 1.8K 



uncertainty. This compares with ou r previous Jupiter estim ate of TjupUer 
and calibration uncertainty of 1.3% (IReadhead et al.ll2004al ). 



3. Data Analysis 
3.1. Improvements to the Analysis Pipeline 

We use a maximum-likelih ood algorithm to c alculate the power spectrum. It is based 



on the framework described in iMyers et al.l (120031 ) (hereafter Paper 4), but includes modi 



fications to deal with the heterogeneous nature of the combined CBI dataset. The pipeline 
consists of two stages: the first stage compresses hundreds of thousands of visibilities into 
~ 10^ "gridded estimators", and calculates noise, CMB signal, and point source covariances 
for those estimators; the second stage calculates the maximum-likelihood spectrum from the 
estimators, as well as ancillary data products. 



3.1.1. Compressing the Data 

As is the case for any compact array, the visibilities in the CBI data have highly cor- 
related signals. To compress the data, we grid the visibilities onto a set of regularly spaced 
estimators in the (M,f)-plane, using the program CBIGRIDR (described in Paper 4). This 
compresses the few hundred thousand visibilities in a typical CBI mosaic into a few thousand 
estimators with essentially no loss of information, under the assumption of Gaussian noise. 
One of the outputs of CBIGRIDR is the matrix R which maps the (Fourier-plane) sky map 
into estimators, 

A = Rt + n, (1) 



^This is the 33 GHz Rayleigh- Jeans brightness temperature of the planet minus the Rayleigh- Jeans 
brightness temperature of the CMB at the same frequency. 
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where t is the Fourier transform of the true sky image, A is the gridded estimator vector, 
and n is the noise component of the gridded estimator (Equation 23 of Paper 4). 

The covariance of the gridded data vector A has signal and noise components 

C = ( A AT) = + + C^'^" + + C''' (2) 

with contributions from noise, CMB, scan-dependent systematic errors, discrete point sources, 
and residual foregrounds (including a "field" of weak, confused point sources too faint for the 
CBI to detect individually) respectively. The f operator denotes the Hermitian conjugate 
(complex conjugate of the matrix or vector transpose). The signal covariance C'^ due to the 
CMB temperature anisotropy signal is 

C^ = RTR^ (3) 

where 

T^(ttt). (4) 

For a Gaussian random CMB temperature field, we expect T to be diagonal if t is represented 
in spherical harmonic set of a^m's, and thus 

{timti'n.') = Ce6{e,£')6{m,m') (5) 

encodes the CMB power spectrum Ci. This remains true in the small-angle approximation 
where 

Tw = {t{ui) t*{ui,)) = Ce 6{ui - u,) (6) 

when t is represe nted in Fourier space and u is the 2-D wavevector, with i = 27i\u\ 
( jWhite et al.lll999l ). Since R contains all the information about the effects of the instrument 
(such as the primary beam) and the scan strategy on the data, to combine heterogeneous 
datasets, we simply grid them separately and then combine their R's. 

In practice, the real and imaginary parts of the complex estimators A are unpacked 
into a real estimator vector d, which is then used for the covariance analysis and calculation 
of the angular power spectrum. 



3.1.2. Calculation of the Power Spectrum 

The log-likelihood of Gaussian data given the model for its signal and noise given in 
Equation [2] is 

log(/:) = -Vc-iA-ilog(|C|). (7) 
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The maximum-likelihood solution is the set of parameters qb which define the theory covari- 
ance Ct {qs) such that the likelihood is maximized. We restrict ourselves to models of Ct 
of the form Ct = YI^bCb- One model of this form is if the Cb are bins in i, in which case 
the qb are the binned power spectrum. 

We calculate the maximum-likelihood CMB power spectrum fror a the es t imato rs using 
MPILIKELY, an MPI implementation of an algorithm described in ISieverd (120031 ) and in 



Appendix \M It requires a single matrix inversion, with no additional operations to 
calculate the gradient and approximate curvature of the likelihood at a trial spectrum. We 
operate on real data, as the correlations for the real and imaginary components of visibilities 
can be different (see iMyers et al.l (120031 )). 



The one algorithmic change we have made in MPILIKELY with the potential to affect 
the output power spectrum is in the numerical treatment of the projection of ground spillover 
and point sources with known positions. To project sources with known positions, but 
unknown fluxes, we form the matrix C'*'''^ = ^i^l i where Si is the expected signal from the 
source, and in the past have added jSC'^^'^ to the noise, for large 0. Similarly, for data 



taken in strips with common ground (i.e., that from ISievers et al.l (120071 )). we calculate the 
matrix C**™" expected from the ground, and have also added ^C'*'^'^"' with a large 7 to the 
covariance. By calculating the spectrum expected from the CMB (see IA.2|) . we flnd that the 
spectrum calculated using large but flnite values for both /? and 7 leads to a recovered CMB 
power spectrum that is biased slightly low in the damping tail due to numerical artifacts in 
matrix inversion - see Figure [H We now take the analytic limit as (3 —>■ 00 (see Appendix 
IA.3|) in MPILIKELY, and explicitly subtract the ground signal as described in Section [21 
with CBIGRIDR accounting for the correlations induced by the subtraction. As can be 
seen in Figure [H the ground subtraction and analytic source projection recovers an unbiased 
CMB power spectrum with slightly reduced errors. 



3.2. Point Sources 

Point sources are the largest astrophysical foreground in the CBI data, and are especially 
important at high-£. All sources with positions that are known reliably from low frequency 
radio observations are removed from our power spectrum analysis; sources below this thresh- 
old require a statistical correction to the power spectrum. An accompanying paper describes 
a cam paign of GBT and OVRO 31 GHz measurements of 2, 125 NVSS sources in the CBI 



fields ( iMason et al.ll2008l ). As a result of this campaign we are able to much more accurately 



determine this correction and characterize its uncertainty. 
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The NRAO VLA Sky Survey fNVSS ICondon et all Il998h is a 1.4 GHz survey of the 
northern sky, taken to be complete down to 5*1.4 = 3.4 mJy (although its nominal detection 
limit is 2.5 mJy). All sources with integrated flux densities above 3.4 mJy in the NVSS 
catalog which lie within one degree of a CBI fleld are projected. Source projection works 
well as long as the relative responses to the source are known for all data. This is true if 
either the source flux is constant for all observations, or all the data are taken in similar 
mod es (same beam, obs erving strategy etc.). Because 30 GHz sources can be quite variable 
(e.g. Iciearv et allboosh and the UV coverage and scan strategy changed between the two 
CBI data sets, we project all sources separately from the scan- and differenced- data. In 
other words, we project two vectors for each source, one corresponding to the source in the 
differenced- data, and one in the scan-data. 

Sources fainter than 3.4 mJy at 1.4 GHz must be accounted for statistic ally. Although 



their sky density is statistically well characterized down to /iJy levels (e.g. iHopkins et al. 



20031 ) from deep observations in a variety of flelds substantially smaller than the CBI flelds, 
these sources are the major systematic uncertainty in the power at high-£ owing to the need 
for a spectral extrapolation from 1.4 GHz, where their counts are well known and where 
they are selected for inclusion in the statistical term, to 31 GHz. Our analysis of GBT and 
OVRO 31 GHz flux density measurements in comparison with the NVSS 1.4 GHz values 
yields an average flux density ratio / = S^i/ Sia = 0.111 ± 0.003; using the low frequency 
counts and the full distribution of / (thus including its intrinsic width as constrained by our 
point source dataset) we determine a point source correction 0.046 ±0.018Jy^/Sr at 31 GHz. 
The probability of an extreme source event giving rise to the high-£ power is also low - the 
highest level of power in a ny of more than 200 realizations of sources consistent with our data 
is less than O.lJy^/sr (see iMason et al.l (120081 ) for the full non-Gaussian distribution). Our 
new source power is nearly a factor of two lower than the earl ier statistical correction used 
in earlier CBI analyses, 0.08 ± 0.04Jy^/sr ( Mason et al.ll2003l ). although consistent within 
uncertainties. The earlier correction was a conservative estimate based on a spectral index 
distribution which was biased against steep-spectrum sources. The current determination 
is based on simulations proceeding from mock CBI observations of the given population of 
residual sources through the full power spectrum pipeline. Amongst other effects this includes 
the effects of Poisson uncertainty in the faint source population. The full distribution of the 
source correction is taken into account when we estimate cosm ological parameters . Further 
details on the point source observations and simulations are in iMason et al.l (120081 ). 



Although our GBT and OVRO data provide strong constraints on the spectral properties 
of the mJy-level extragalactic sources that dominate the CBI statistical correction, it is 
possible that at fainter flux densities the sources have different characteristics at 30 GHz. The 
majority of the fainter 1.4 GHz sources are expected to have steep radio spectra (IRichards 
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2000l : ICondon et al.lll998l ). however, were there to be a substantial enhancement of sources 



with strong ly inverted spectra between 1.4 and 31 GHz they could contribute appreciably to 
the excess. iMason et al.l ( 20081 ) considers this scenario in detail, finding that for moderately 
inverted spectra (a ~ 0.2) they would need to constitute 40% of the sub-mJy population in 
order to fully explain the CBI excess. Were the sources to have strongly inverted spectra 
(a ~ 0.8), 2% of the population would be required. In contrast, in the GBT+OVRO surveys, 
the most steeply inverted spectrum source had a = 0.49 and < 0.1% of sources had a > 0.3. 



3.3. Data Tests 



Th e data presented here hav e been discussed p r eviousl y in lMason et al.l (120031) .iPearson et al 



(]2003h . iReadhead et all (l2004ar ). iReadhead et al.l fl2004bf ). and ISievers et all (120071 ). where 
extensive data integrity tests were described. A number of further tests have been carried 
out on the data. Key results from this exercise are as follows: 



3. 



4. 



Dish Pointing Errors: By analyzing beam maps on bright sources we have deter- 
mined the individual antenna primary-beam pointing errors to be ~ 3.5' RMS in a 
single direction. We have simulated the impact of this, folding in the correlations in- 
duced by deck rotations throughout our observing strategy along with our typical 15" 
RMS pointing errors, and find that the resulting bias in the power spectrum is < 9 yuK^ 
in the highest i bin. 

Spectral Index of Projected Sources: We have introduced systematic errors of 
6a = ±1 in the spectral index used to project sources out of the data and find that 
this has an effect < 6 /iK^ in the highest £ bin. 

Primary Beam: To test the sensitivity of our power spectrum to the precise beam 
shape used in the analysis, we ran simulatio ns of CBI observatio ns using our best-fit 
physical model of the CBI beam (described in lPearson et al.ll2003l ). but analyzed them 
with the best-fit Gaussian beam. This results in < 10 fiK^ error at any i. 

Noise Fitting: The way that we calculate the thermal noise for individual visibilities 
in our dataset results in a known ~ 2 — 6% underestimate of the thermal noise variance. 
This underestimate depends on the observing strategy that was used. Sirnulatiori s and 
analyti c calcu lations of this effect have been presented in iMason et al.l (120031 ) and 



Sieverd (120031 ). To further constrain the thermal noise power spectrum, which can be 
a limiting factor especially at the highest i CBI measures, we have implemented a 
noise estimator after the gridding step in the analysis pipeline. This estimator splits 
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the 10 CBI frequency channels into two groups [e.g., low-frequency/high- frequency, 
or even/odd channels) and fits for a thermal noise multiplier and the CMB power 
spectrum simultaneously. For the polarization observations this exercise yields noise 
spectrum multipliers of 1.0195 + / — 0.009 (low/high) and 1.015 + / — 0.009 (even/odd), 
in comparison with our previous best estimate of 1.0175. For the total intensity obser- 
vations we find 1.071 + /0.008 and 1.071 + / — 0.008, in comparison with our previous 
best estimate of 1.057. We adopt noise variance multipliers of 1.017 and 1.071 for the 
polarization and total intensity observations, respectively, and an uncertainty of 0.9%. 

5. Bright Sources: The large (~ 143 deg^) area covered by CBI makes it impossible 
to completely avoid bright radio sources: there are two 5*30 > 300 mJy sources in the 
fields. For the observed distribution of sources in flux density (A^(> S) ~ S~^), the 
brightest sources (or constant fractional residuals to them) will dominate the map 
variance. To measure the possible effect of these sources we reanalyzed the data after 
removing all individual CBI pointings where a point source with an apparent flux 
density 6*30 > SOmJy was evident, which removed 33 out of 259 pointings. This 
reduces the power spectrum by less than 15 /iK^ averaged over all i and shows no 
characteristic trend of increasing to high-£, which would be expected of residual source 
contamination. 

6. Source Structure: We project pure point source templates for the NVSS sources. 
The NVSS resolves ~ 20 percent of the sources it detects. If the sources are resolved 
by the CBI, then we would expect leakage into the power spectrum. Using Montag^ 
we mosaic the NVSS 4 deg x 4 deg maps to cover each of our CBI fields, and, zeroing 
out all pixels below 3 times the NVSS RMS noise level, simulate CBI observations of 
those maps. These observations are run through the full CBI data reduction pipeline, 
projecting out sources using the NVSS catalogs above 3.4 mJy. We find that our 
projection method removes >99.9% of the source power, with the residual power a 
factor of 20 times lower than the observed CBI high-£ signal. This is a conservative 
estimate of the signal from source extendedness, since the high-frequency emission 
from AGNs tends to preferentially come from compact cores. In addition to testing 
the impact of source size, this is a powerful test of our data pipeline and the quality 
of the NVSS catalog. 



^ lit tp : / / montage . ipac . caltech . edu 
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4. Results & Interpretation 

4.1. Power Spectrum 

The final CBl total intensity power spectrum is shown in Figure [H together with the 
WMAP 5-year, ACBAR, and QUaD power spectra, and best-fit tilted ACDM power spec- 
trum model. The spectrum is given in Table [2l and full window functions and bin-bin 
correlations are available on-linejj 

A comparison of these results with the previous CBI power spectrum is shown in Fig- 
ure O The marked improvement is due to several factors: the inclusion of ~ 50% more data, 
use of the GET 30 GHz observations to reduce the uncertainty due to point sources, and 
the algorithmic improvements described in § 13.11 

At £ > 1800 there remains a clear excess of power ove r the expected intrins ic CMB 
anisotropy. This i s shown in Fi gure HI along with A CBAR (IReichardt et al.ll2008l ). QUaD 
( Prvke et al. l20o"8), and BIM A ( Dawson et al. 2006) measur e ments . To quantify this we 
use the iKomatsu &: SeljakI (120021 ) analytic and the iBond et al.l (120051 ) SPH simulation-based 
predictions of the SZ angular power spectrum. (See Section 14.31 for more details.) An 
amplitude for the SZ spectrum is then fit in conjunction with an amplitude for the intrinsic 
anisotropy spectrum, using a canonical tilted-ACDM mode£| as a shape. The templates are 
calculated with cxs = 0.8 and Qbh = 0.0321. This results in a 3. la detection of power in 
excess of the intrinsic anisotropy (a best-fit scaling of 3.5 ± 1.3 of the nominal KS template). 
The SPH template has a best-fit scaling of 5.4 ±1.8 and a 3.4cr detection significance. The 
detection significances are calculated using \/—26 log (£), where the likelihood difference is 
between the best-fit CMB-only shaped model and the best-fit CMB+template model. We 
show the total power spectrum including the contribution of the best-fit SZ spectrum by the 
dashed lines in Figure HI Of note is the marked change in the magnitude of the SZ power 
spectrum between 30 and 150 GHz, and the fairly broad contribution of clusters down to as 
low as £ ~ 1000. 



Since the secondary SZ anisotropy will be highly non-Gaussian, uncommon structures in 
the ~ 4 deg^ of "deep" CBI pointings could bias the power spectrum (although the statistical 
weight of the wide, shallow area is about 2.5 times of that of the deep data). We re-ran the 
power spectrum extraction both excluding the deep field data from the analysis, and using 



^http:/ /www. astro. caltech.cdu/'^tjp/CBI/data/index. html 



''Flat, nbh'^ = 0.0223, flcoM = 0.108, r = 0.087, Us = 0.96, and including the effects of gravitational 
lensing. No SZ contribution is included, as that is handled separately. 
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only the deep data. The results are shown in Figure [5l While the error bars are increased by 
~ 25%, the amplitude of the power spectrum at ^ > 1000 is not reduced, indicating that the 
small-scale excess power is not a peculiar property of the deep fields observed by CBI. An 
assessment of the cosmic variance in the high-£ CBI power spectrum, under the assumption 
that the dominant signal over intrinsic anisotropy is due to SZ clusters, is presented in § 14.31 

Power spectra from the individual CBI fields are shown in Figure [61 To increase the 
signal-to-noise ratio a coarser binning has been used. Each of the CBI mosaics shows an 
excess of at least 0.7a , with levels of ( 2.7±2.6, 7.6±2.9, 2.6±2.3, and 1.6±2.2 ) above the 
CMB using the as = 0.8 KS template for the (02, 08, 14, 20)-hour mosaics. The individual 
fields are all consistent with the same excess level (reduced = 1-01 per dof), with the 
most discrepant field 1.6a away from the mean. 

The CBI data by themselves cannot determine the source of the excess. In particular, 
they cannot distinguish an SZ foreground from point sources (if one fits for a point source 
and an SZ amplitude simultaneously, the power spectra are sufficiently similar that the error 
bars more than triple). If instead of an SZ template, we model the excess with a point 
source-like template, we find its total amplitude is 0.180 ± 0.052Jy^/sr above the expected 
contribution from unresolved sources of 0.046JyVsr. To high accuracy, the best-fit SZ excess 
level is linearly dependent on the input source level. For the KS template, the best-fit SZ 
amplitude is ^^2=4. 59-24. Txgj^o where g^z is in units of the predicted SZ signal for as = 0.8 
and Qiso is the faint source contribution in Jy^/sr. One can use this relation to ask how 
different mean source levels would affect the excess: a value of 0.186 for g^so would make the 
high-£ excess disappear, a value of 0.146 would make the best-fit CBI value be equal to the 
predicted KS level, and a value of 0.096 would make the best-fit CBI value be within 1-a of 
the predicted KS level. 



4.2. Diffuse Foregrounds 

While on small scales, point sources are the largest 30 GHz foreground, on larger 
scales (^ 1 degree), diffuse foregrounds dominate. At frequencies of ~ 30 GHz, the ma- 
jor known diffuse Galactic foregrounds are synchrotron radiation and free-free emission. 
Vibrational (thermal) dust emission is negligible at frequencies below ~ 50 GHz. How- 
ever, there is considerable evidence for an additional component which is closely corre- 



10-50 GHz ( 


Leitc 


1 et al. 


1997 


Bandav et al. 


2003; 


de Oliveira-Costa e 


t al. 


2004 


Finkbeinei 


2004 




Davies et al. 


20061: 


Bonaldi et al. 


2007; 


Hildebrandt et al. 


2007 




Dobler & Finkbeinei 


2008a 


b). Both the COBE-DMR and WMAP datasets have shown that the high latitude 
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sky is dominated by such a foreground closely linked with dust emission, usually traced in 
the FIR (A ~ 100 fim). The most popular candidate for this an omalous component is th e 
so-called "spinning dust" emission (iDraine fc Lazarian Il998al jbl: lAli-Hai'moud et al.l 120081 ) . 
but the situation is far from clear. 

Whether or not the dominant diffuse emission is spinning dust or flat-spectrum syn- 
chrotron is not of great importance in the context of this paper. However, we need to 
quantify the level at which diffuse foregrounds are present in our data. We use the fact that 
the FIR maps are a good tracer of the foreground morphology at these frequencies. In the 
WMAP data at K and Ka-bands (23 and 33 GHz, respectively), there is a strong correlation 
with the 100 /xm lSchlegel et al.l (119981 ) map. If the emissivity of the dust emission is roughly 
constant for a given region, cross-correlation between the FIR map and CBI data provides 
a very sensitive method to detect such emission. Furthermore, IRAS data have adequate 
resolution to cover the angular scales measured by CBI. 



We simulated CBI data based on the ISchlegel et al.l (119981 ) 100 /im map as our fore 



ground "reference" map. The 100 /im maps were converted to CBI visibilities using the 
MOCKCBI software based on the real observed visibility data sets. To convert to the approx- 



(Bandav et al. 


2003; 


Davies et al. 


2006) 



ity can vary by a factor of a few over the sky, it provides an initial guess for the amplitude 
of the signal i.e., we expect a cross-correlation coefficient of ~ 1 based on previous data on 
la rger angular scales. We also repeated the analysis using Ha data from the compilation 



of iFinkbeineii (120031 ). Ha is known to b e good tracer of free-f ree emission at high Galactic 



latitudes where dust extinction is small (jPickinson et al.ll2003l ). We scaled the Ha template 
by 5.83 fiK per Rayle igh, which is the expected value at 31GHz, assuming Tg ~ 8000 K 



(IDickinson et al.l 120031 1. 



The simulated foreground visibilities were fitted to th e CBI data us i ng the template- 
fitting method of Appendix IA.5I For the dust template of ISchlegel et al.l (119981 ) . we found 
a combined correlation coefficient of 1.18 ± 0.53 (2.2o") for the entire dataset. This suggests 



that a very small level of contamination from Galactic emis sion might exist in t 



le CBI data. 



20031: 



Davies et al. 



at a l evel similar to those observed at larger angular scales ( iBanday et al 

2006 ). The fluctuation power Ci in the FIR dust scale as Ce oc ( Gautier et al. 19921 ) 
thus the power will be mostly at large angular scales. For the Ha template we flnd a 
combined correlation coefficient of 0.78 ± 0.58. For both foreground templates, we find a 
negligible impact on the power spectrum, with the high-£ excess changing by less than 1%. 
The emissivity factor is close to what is expected at high latitudes and indicates that the 
CBI fields are relatively low in foregrounds on angular scales < 1°. This is also supported by 
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the field- by-field data splits ( §4.ip : all the independent fields give consistent power spectra. 



4.3. Cosmological Parameters 



We can place constraints on the standard parameters of tilted ACDM cosmology by 
fitting a range of model s pectra to the obs erved CMB (and other) data. Model spectra 
were generated by CAMB (ILewis et al.ll2000l ) and the parameter const raints determined by a 
modified version of the Monte-Carlo Markov Chain code, COSMOMC flhewis k Bridyl2002[ ). 
For all parameter runs, we take into account the effects of CMB lensing, and assume a prior 
that the universe is spatially flat. To capture the non- Ga ussian natur e of th e CBI's power 
spectrum, we use the offset-lognormal approximation of iBond et al.l (120001 ) in parameter 
analysis. To convert model spectra to predicted CBI bandpowers, we evaluate the window 
functions with a spacing of 51 = 20, and use cubic Hermitian polynomial interpolation 
between the measured points. 



We extend the treatmen t of the SZ foreground relative to previous analyses (IBond et al. 



20051 : iReadhead et al.ll2004al ). We assume the angular powe r spectrum from clusters scales 
in a simple analytic fashion: oc (j| [QbhY (IBond et al.ll2005l : iKomatsu fc Seliakll2002f ). We 
have explored two sets of SZ templates: 



1. Semi-analytic templates from lKomatsu Sz Seljak J2002h (the KSSZ template), updated 
to include the effect of ag on the shape of the SZ power spectrum (E. Komatsu, private 
communication). The dependence of the shape on erg is found to be minimal. 



2. Power spectra from tree-SPH simulations described in iBond et al.l (120051 ) (the SPH 
template). 



The SPH template has relatively less power than KSSZ for flxed cosmological parameters, 
and rises more quickly with i. Unless stated otherwise the KSSZ template is used in all 
parameter analysis. Also unless otherwise indicated, all runs use the full i range of data and 
explicitly model the SZ foreground. Based on a comparison of these theoretical templates, 
we estimate that for a given background cosmology there is a factor of ~ 2 uncertainty in 
the SZ power spectrum. This corresponds to a ~ 10% systematic uncertainty in ag- Because 
of this we do not link the SZ spectra to the background cosmology, e.g., through Qf, or ag. 
Rather we use an independent parameter crf^ to describe the amplitude of the SZ power 
spectrum and to index the family of shaped SZ power spectra. 



Table [3] sh ows the marginalized individual param eter results for an analysis including 
WMAP 5-year jNolta et allboosl : lOunklev et aPboOsh and CBI data, using both SZ tem- 
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plates. This includes a marginalization over the uncertainty in the residual point source 
power spectrum using the distribution determined from the simulations in § 13.21 Table H] 
contains t he same parameters, b ut with t he addition of more CMB d ata to CBI and WMAP : 
ACBAR jReich^rdjejallboOsh. BIMA JPawson et al.l 120061 ) . VSA Jpickinson et al.l[20oi ). 
Boomerang jMontrov et all l20od : I Jones et all l2006h . QUAD jAde et al.l l2008l. adopting 
their "pipeline 1" spectrum), and both CBI and DASI polarization data (jsievers et al. 2007 : 
Leitch et al.ll2004j ). We henceforth refer to this data combination as CMBall. The choice of 
SZ templat e has a n impact on the value of erg inferred from the primary fluctuations - the 
Bond et al.l (120051) template gives ag values that are systematically higher by about 0.015 



than the iKomatsu &: Seljakl (|2002| ) template. This is because the KSSZ template is flatter 
than the SPH template, so for fixed observed power at ^ ~ 2000, KSSZ will remove more 
power il l the region of the th ird peak, dropping ag. Our values of o"g are also lower than 
those in iDunkley et al.l (120081 ) where the level of the KSSZ template is capped at twice the 
level predicted for a ug = 0.8 universe. We allow the amplitude to float freely, and the CBI 
data set the level to be around 2.5 times the ag = 0.8 prediction, again resulting in more 
power being removed from the third peak and a lower primary ag value. 

The values of cxg inferred from the high-£ power spectrum are af^ = 0.910 ± 0.064 for 
CBI+WMAP5 and af^ = 0.922 ± 0.047 for CMBall. As a point of comparison, using the 
SPH template, we find af^ = 1.015 ± 0.060 for CBI+WMAP5 and af^ = 0.977 ± 0.049 
for CMBall. CMBall has a (slightly) higher inferred ap than CBI+WMAP5 for the KSSZ 
template, whereas it has a lower one for the SPH template. This is due to the BIMA 
data falling above the CBI prediction for the relatively flat KSSZ template, while they fall 
below the CBI prediction for the SPH template. We find that CBI, ACBAR, and BIMA are 
consistent with each other if the high-£ excess is due to SZ - see Figure [7] for full MCMC 
chains comparing the three experiments using both SZ templates. For comparison with other 
experiments, we give the best-fit CBI excess level (from Table E]) for the two templates at 
selected £ and frequencies. The results are summarized in Table O 

When we extend the ACDM model to include tensor modes, we find no detection of 
them: r < 0.32 (95%) for CMBall, where r is the tensor-to-scalar power ratio at the pivot 
wavenumber 0.05 h Mpc~^. 

We ha ve explored the poten tial impact of cosmic strings in the CMB using a string 
template of iPogosian et al.l (120081 ). We treat the string contribution as a template of fixed 
shape added to the power spectrum, and allow its overall amplitude qstring to float, just as we 



have for the SZ templates with overall parameter qsz- We re-emphasize the iPogosian et al. 



( 120081 ) caution that the power spectrum from strings is not uniquely determined and its 
shape depends on the details of the string properties. While no combination of CMB data 
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detects cosmic strings, the addition of high-£ CMB data adds a significant constraint to the 
maximum allowed string amplitude: the 95% string upper limit from CMBall is 65% of the 
limit from WMAP5 alone (see Fig. [H]). The high-£ data constrain strings because the power 
spectrum from them generically falls off more slowly with I than that fro m the adiabatic 



fluctu ations since the string fluctuations are not subject to Silk damping (jPogosian et al. 
20081 ) ■ The amplitude Qstring oc {Giif can be expressed in terms of the string tension ji 



times Newton's constant G. The limits on G/i (which scales like the square root of the 
power spectrum amplitude) are G/i < 3.4 x 10^^ for WMAP5 only, and G/i < 2.8 x 10^'' 
for CMBall. With WMAP only, the string amplitude is highly degenerate with Us {e.g., 



Battye et al.l (120061 ) . see Fig. [8]), and there is no longer a "detection" of < 1, rather 
we obtain Hs = 0.972 ± 0.018 with a 95% upper limit of 1.007. The addition of the high-£ 
CMB breaks this degeneracy, and remains less than one with high significance: we obtain 
ris = 0.961 ± 0.015 with a 95% upper limit of 0.990. 

We find that the addition of the CBI data have a strong impact on the running of 
the scalar spectral index, -^^^ (evaluated at a pivot wavenumber of 0.05 h Mpc~^). For 



the KSSZ template, we find = (-0.041 ± 0.031, -0.048 ± 0.028, -0.066 ± 0.022) for ( 
WMAP5, WMAP5+CBI, CMBall), respectively. For the SPH SZ template, we find = 
(-0.039 ± 0.030, -0.042 ± 0.027, -0.059 ± 0.022) for ( WMAP5, WMAP5+CBI, CMBall), 
respectively. Thus, with CMBall, the running of the spectral index is a 2.76-0" "detection" 
for the SPH template, and 3.03-0" one for the KS template. The detection is substantially 
driven by the high values obtained for the SZ signal. When the allowed SZ level is capped 
at twice the nominal KS value, we find that limits on for CMBall are —0.048 ± 0.021 

' a log fe 

and —0.043 ± 0.021 for the KS and SPH templates, around 2-o" detections (see Fig. [H]). We 
stress that the inferred value of -^^^ depends not only on the amplitude of the SZ signal, 
but somewhat on its shape as well, hence the difference in between the two templates. 

Since clusters are compact sources, the distribution from field to field will not be Gaus- 
sian. We have assessed the effect on the sample variance i n the SZ clust er power spectrum 



for the CBI coverage region using the simulated maps of IWhitd (120031 ). These consist of 
10 maps of the thermal SZ effect, each 10° x 10° in size. These are calculated using large 
N-body simulations in a ACDM cosmology with ag = 1. The pressure profiles of the clusters 
were made assuming the gas density foll ows the dark matter d ensity, and the temperature 
is isothermal and proportional to M^^f^ ( Schulz fc White 2003). The SZ maps are line-of- 



sight integrations of these 3D pressure configurations. Fake observations of these fields are 
constructed using the real CBI uv coverage and per-visibility noise levels and they are run 
through our power spectrum extraction. The thermal SZ power spectrum at £ > 1800 emerg- 
ing from this analysis has a mean level of ~ 200 /xK^. When only the CBI deep fields are used 
we find an RMS in the SZ power spectrum of ~ 100 /iK^, that is a fractional scatter of about 
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50% in the (noiseless) spectra. When the full CBI sky coverage is probed instead of just 
the deep fields, we find a fractional scatter in the SZ spectrum to be 21%. The 21% scatter 
could be driven by non-Gaussian errors in the deep fields. To test this, we excluded the deep 
fields from the CBI mosaics and repeated the analysis. We found that the fractional scatter 
dropped to 19%, hence the deep fields do not drive up the variance of the total result by a 
substantial amount. We expect that the non-Gaussian component will be quite sensitive to 
CTg since the number of clusters drops dramatically as erg is dropped, increasing the Poisson 
fluctuations. Thus the sample variance in the SZ spectrum depends heavily on the value of 
(78, and since the sample variance errors are highly correlated between multipole bins, these 
results have not been fed into further analysis, but should be taken as indicative. 



5. Summciry 

We have presented the total intensity power spectrum resulting from five years of ded- 
icated CBI observations and campaigns of point source foreground observations with the 
OVRO 40-m and the GBT. The OVRO and GBT data allow us to greatly improve our esti- 
mate of the power from faint 30 GHz radio sources. On its own, the CBI cannot distinguish 
between power spectra from point sources and from other sources such as the SZ effect from 
galaxy clusters, and so these supporting observations are essential. Our data support the 
existence of excess power above the primary CMB on small angular scales at ~ 3a, and are 
the most sensitive constraints to date on the statistical SZ cluster foreground. In extensive 
testing, we find no evidence of any instrumental systematic effect capable of giving rise to the 
excess. By running the NVSS maps (with all pixels below three times the NVSS RMS noise 
zeroed out) through the CBI pipehne, we find that our source treatment rejects >99.9% of 
the power in known sources. This test confirms the quality of the NVSS catalogs, the vahdity 
of the projection of known sources, and the insensitivity of the CBI to extended sources. If 
the excess is due to the SZ effect, we find that other data, notably ACBAR and BIMA, 
are consistent with the excess seen by CBI. For two different SZ spectral templates, we find 
that the inferred (jg from the excess is marginally inconsistent with those derived from the 
primary fluctuations, 0.922 ± 0.047 and 0.988 ± 0.049 vs. 0.769 ± 0.031 and 0.784 ± 0.030. 
To determine definitively the imphcations of the observed small-scale excess power further 
theoretical work is needed. We find that high-£ data break the degeneracy between the tilt of 
the spectral index Ug and potential contributions from cosmic strings. Prom running n-body 
simulations through the CBI pipehne, we find that the scatter in the expected level of the 
signal from clusters due to their non-Gaussian nature is about 20%. 
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Fig. 1. — Effects of projections on the expected spectrum from tlie CBI strips. Tlie red 
triangles sliow tlie expected spectrum from just the TT part of the CBI polarization data 
with no removal of ground or sources. The green squares show the expected spectrum when 
the scan mean (which contains the ground signal) is projected out. The red X's show the 
spectrum when only the sources are projected. The green X's show the spectrum when both 
sources and ground are projected - there is a significant bias downwards in the spectrum. 
The grey circles show the the same spectrum when the scan mean is subtracted rather than 
projected, and the sources are projected using the numerically stable method of Appendix Rl 
The black curve is the input CMB spectrum. This method keeps the covariance matrix better 
conditioned, so it is less susceptible to roundoff errors in the inversion. These spectra are 
fit to the theoretical CMB+noise signal matrix using the techniques described in Appendix 
\A.2\ which is equival ent to averaging oyer all possible noise and signal simulations. The CBI 
differenced data from lReadhead et al.l (l2004al ) were already effectively scan-mean subtracted, 
and so are not included here. 
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Fig. 2. — CBI total intensity power spectrum. The blue points are the CBI power spectrum 
in this wo rk, given iri text form in Table [2l The salmon points are the W MAP 5-year 
spect rum (INolta et al.l 120081 ). The green points are the QUaD 2-year spectrum (jPrvke et al. 



20081 ). The burnt sienna points are the ACBAR 150 GHz spectrum flReichardt et al. 



20081). 



The full CBI spectrum, including bin correlations and window functions, is available online. 
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Fig. 3. — A comparison of the CBI 5-year power spectru m, including results from the GBT 
30 GHz survey, with the two-year CBI power spectrum of iReadhead et al.l (l2004al ). Note the 
changed finer b inning at high i f or the power spectrum from the 5-year data. The window 
function of the iReadhead et al.l (l2004al ) highest-^ bin extends from i ~ 2000 to £ ~ 3500. 
Note that the error bar on that very b ig bin is about the same as for the last of the finer bins, 
in spite of the large IReadhead et al.l (l2004al ) bin being broken up into three distinct bins in 
this work. The main reasons for the improvement in the spectrum is the factor of two more 
data, and the development of analysis techniques that allowed us to combine these disparate 
datasets. In the dampin g tail, the new spectrum is about a factor of two improvement over 
R,eadhead et^ fl2004al ). 



-22- 




Fig. 4. — The CBI, ACBAR, QUaD, and BIMA power spectra at small angular scales 
are contrasted. The solid black line shows the tilted ACDM model from Section 14.11 for 
the CMB primary anisotropi cs. The dashed line s inclu de the contribution of secondary SZ 
anisotropy using the model of lKomatsu &: SeljakI (l2002l ) with the best-fit template scaling of 
3.5 (in bandpower) that we have determined using only the CBI data. The SZ plus primary 
anisotropy power combination at 30 GHz is the blue-dashed line. Note th at, apart from 
fitting the CBI power spectrum, it passes through the BIMA point at £ ~ 5300 (jPawson et al. 



20061 ). We have also forecast the level for SZ plus primary anisotropics at 150 GHz (red 



dashed line) and 1 00 GHz (orange dot-d ash line). These are compar ed with the power 
spectra of ACBAR jReichardt et ahlboosh and QUaD JPrvke et al.lboosh . 
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Fig. 5. — A comparison of the power spectrum obtained from all CBI TT data with that 
obtained when the deep fields are excluded, and with that obtained when only the deep 
fields are used. Although the error for the no-deeps at £ ~ 2500 is larger, and consistent 
with no excess at the ~ 1 — a level, the overall mean amplitude is about the same as for the 
deep-only case. The all-data and no-deeps spectra are at the same i, but have been offset 
for clarity. 
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Fig. 6. — Power spectra from individual CBI fields show how the bandpowers fluctuate from 
field to field. The individual fields are all consistent with each other, and each sees power 
above the CMB. The spectra are staggered in £ for clarity in plotting. 
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Table 1: CBI Fields 
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Fig. 7. — 1- and 2-cr likelihood contours of the SZ amplitude and the baryon density, as 
determined from MCMC chains for CBI, ACBAR, and BIMA, as indicated. We ran chains 
for WMAP5 plus, in turn, CBI, ACBAR, and BIMA, fitting an SZ excess template to each 
case. The left panel shows the results when using the KSSZ template, and the right panel 
shows the same using the SPH SZ template, marginalizing over the other parameters. For 
both templates, the \-a regions of the three experiments are all in excellent agreement if the 
excess is due to SZ. 
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Table 2: CBI Power Spectrum. Power spectrum from the total CBI dataset. Columns are 1) 
bin, 2) power spectrum, 3) Gaussian error on the power spectrum, 4) thermal noise power 
in the bin, 5) the contribution from unresolved source in the bin (assuming 0.046 JyVsr), 6) 
lower £ limit, and 7) upper ^ limit. 





Q(/xK2) 


Error 


Cjioisei^l^^ ) 




^rmn 




1 


5695.372 


753.175 


1767.465 


2.4 


0.00 


350.00 


2 


1260.879 


260.277 


218.186 


1.2 


350.00 


450.00 


3 


2987.370 


408.622 


300.134 


3.5 


450.00 


540.00 


4 


2449.122 


376.215 


235.321 


4.7 


540.00 


620.00 
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1844.757 
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311.309 


3.5 


620.00 


700.00 
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4.7 


700.00 
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860.00 
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940.00 


9 
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221.538 


600.919 


9.9 
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10 


979.649 


220.328 


555.848 


9.5 
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11 


1166.825 
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607.624 


11.9 
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1192.00 


12 


1041.110 


166.131 


669.178 


14.2 


1192.00 


1302.00 


13 


727.473 


143.316 


706.893 


13.7 


1302.00 


1425.00 


14 


803.088 


147.134 


1218.974 


19.7 


1425.00 


1600.00 


15 


311.402 


126.647 


1240.392 


27.4 


1600.00 


1800.00 


16 


379.102 


122.922 


1316.752 


29.9 


1800.00 


2050.00 


17 


260.743 


132.234 


1666.893 


35.4 


2050.00 


2350.00 


18 


386.738 


117.427 


3000.407 


66.1 


2350.00 


3900.00 
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Tablc 3: Cosmological Parameters from CBI+WMAP5. These are the standard basic 6 
parameters of the tilted ACDM model and the SZ amplitude qsz, plus parameters derived 
from them for CBI and WMAP5. Note that all of the parameters except for as, and the as 
inferred from the SZ amplitude erf ^ oc q^^ are relatively insensitive to the two SZ templates 
used. All parameters above the line are independent variables, and those below are derived 
from the independent parameters. 



Parameter 


CBI+WMAP5+KS sz 


CBI±WMAP5±Bond sz 




0.02291 ±0.00061 


0.02271 ± 0.00060 




0.1069 ±0.0064 


0.1081 ±0.0063 


e 


1.0406 ± 0.0030 


1.0404 ±0.0030 


T 


0.087 ±0.018 


0.086 ±0.018 


qsz 


2.52 ±0.96 


5.3 ± 1.8 


ris 


0.960 ±0.015 


0.963 ±0.015 


log[10^^A,] 


3.039 ± 0.044 


3.052 ± 0.042 




0.756 ± 0.030 


0.750 ± 0.030 


Age/GYr 


13.65 ±0.14 


13.68 ±0.14 




0.244 ±0.030 


0.250 ±0.030 


CTS 


0.770 ± 0.038 


0.781 ±0.036 


Zre 


10.8 ± 1.5 


10.9 ± 1.5 


SZ 


0.910 ±0.064 


1.015 ±0.060 


Ho 


73.2 ± 2.8 


72.6 ±2.8 
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Table 4: Cosmological Parameters Determined from the Combined CMB Datasets. These are 
the standard basic 6 parameters of the tilted ACDM model and the SZ amplitude qsz^ plus 
parameters derived from them for combined CMB datasets (see text for the list). Columns 



I—' O T* O /~\"t~ 

jT ai ametei 


vv ivi/\j:^o-rvyivii3aii-|-ivo oZj 


VV iVl/\jr 0-rvjiVii3ail-rojr rl oZj 




u.uzzoy It u.uuuo( 


U.UZZ04 It u.uuuoo 




U.iU(0 It U.UUOo 


U.iUoo It U.UU04 


e 


1.0419 ± 0.0025 


1.0417 ±0.0025 


T 


0.085 ±0.018 


0.085 ±0.017 


OiSZ 


2.65 ±0.81 


4.3 ± 1.3 


ris 


0.956 ±0.014 


0.960 ±0.014 




3.033 ±0.039 


3.051 ±0.038 




0.758 ±0.026 


0.749 ± 0.027 


Age/GYr 


13.62 ±0.13 


13.66 ±0.13 




0.242 ±0.026 


0.251 ±0.027 


0-8 


0.769 ±0.031 


0.784 ±0.030 


^re 


10.7 ± 1.5 


10.8 ± 1.5 


SZ 
'^8 


0.922 ±0.047 


0.988 ± 0.049 




73.5 ±2.6 


72.6 ±2.6 



Table 5: CBI Predictions for Other Experiments. For convenience in comparing the CBI 
results with other experiments if the excess power is due to SZ, we present the signals 
predicted by the two SZ templates at a variety of frequencies and angular scales. All values 
are in /xK^. 





30 GHz KS 


100 GH 


z KS 


150 GHz KS 


30 GHz SPH 


100 GHz SPH 


150 GHz SPH 


500 


48.6 


± 


18.5 


29.0 


± 


11.0 


11.6 


± 


4.4 


24.3 ± 8.3 


14.5 ± 4.9 


5.8 ± 2.0 


1000 


73.3 


± 


27.9 


43.6 


± 


16.6 


17.4 


± 


6.6 


52.5 ± 17.8 


31.3 ± 10.6 


12.5 ± 4.2 


1500 


87.8 


± 


33.4 


52.3 


± 


19.9 


20.9 


± 


8.0 


94.2 ± 32.0 


56.1 ± 19.1 


22.4 ± 7.6 


2000 


94.2 


± 


35.9 


56.1 


± 


21.4 


22.4 


± 


8.5 


138.2 ± 46.9 


82.3 ± 28.0 


32.9 ± 11.2 


2500 


96.2 


± 


36.6 


57.3 


± 


21.8 


22.9 


± 


8.7 


151.8 ± 51.5 


90.4 ± 30.7 


36.1 ± 12.3 


3000 


96.2 


± 


36.7 


57.3 


± 


21.8 


22.9 


± 


8.7 


189.3 ± 64.3 


112.8 ± 38.3 


45.0 ± 15.3 


4000 


92.4 


± 


35.2 


55.1 


± 


21.0 


22.0 


± 


8.4 


232.9 ± 79.1 


138.8 ± 47.1 


55.4 ± 18.8 


5000 


86.2 


± 


32.8 


51.4 


± 


19.6 


20.5 


± 


7.8 


280.4 ± 95.2 


167.0 ± 56.7 


66.7 ± 22.7 


6000 


80.3 


± 


30.6 


47.9 


± 


18.2 


19.1 


± 


7.3 


294.0 ± 99.9 


175.2 ± 59.5 


70.0 ± 23.8 


8000 


67.4 


± 


25.7 


40.1 


± 


15.3 


16.0 


± 


6.1 


316.9 ± 107.6 


188.8 ± 64.1 


75.4 ± 25.6 
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WMAP5 

WMAP5+CMBall 



1.02 



Fig. 8. — 1- and 2-ct likelihood contours of cosmic string amplitude and n .^, marginalized 
over o ther parameters. The string amplitude is relative to the Q-template of iPogosian et al. 
(120081 ) . which is normalized to a string tension of Gfi = 1.1 x 10~^. With only WMAP5 
data, Us is partially degenerate with qstrmg, and ra^ < 1 is no longer significant at a 2-cr level. 
The addition of the high-^ CMB data breaks this degeneracy, and < 1 at the 2-a level 
holds even with the addition of cosmic strings to the parameter analysis. 
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Fig. 9. — 1- and 2-cr likelihood contours of the running of the spectral index rirun = 
■^^^{0.05h Mpc~^), and SZ amplitude qsz, marginahzed over other parameters. A high 
level of the SZ template pulls power out of the primary CMB fluctuations in the region of 
the third peak. This drives more n egative. The a r tificia l cropping of qsz at a value 

of 2.0 as used in the WMAP5 analysis of iDunkley et al.l (120081 ) clearly distorts the picture 
relative to a freely floating qsz- 
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A. Maximum Likelihood Fitting 



We use the program MPILIKELY to measure the maximum-likelihood power spectrum 
and related quantities from the gridded data and its no ise and s ignal matrices. MPILIKELY 
is an MPI implementation of the algorithm described in lSieverd (120031 ) . We briefly summarize 
the algorithm here, as well as describe additional features of MPILIKELY. 



A.l. Fast Curvature and Gradient Calculation 



The likelihood of correlated Gaussian data is: 



log(/:) = -idtC-id-ilog(|C|) (Al) 



where d is the data, and C = {didj) is the covariance matrix( e.g., iBond et al.l (119981 )). 
The correlation matrix in general will depend on both the noise and the signal in the data. 
The maximum-likelihood solution is then the set of parameters on which C depends that 
maximizes the likelihood. It is often the case that C depends linearly on its parameters. 
This is true if we parameterize the CMB power spectrum by bands in £, in which case the 
theory matrix of Equation [2] takes the form 



Ct = Y,^bCb, (A2) 



in terms of CMB signal matrices Cb with associated bandpowers g^- The standard technique 
for maximizing the likelihood is to calculate the gradient and curvature of the log-likelihood 



This preprint was prepared with the AAS IMJjjX macros v5.2. 



- 35 - 



and take a mult i- dimensional Newton's method step, iterating until convergence. The gra- 
dient of the log likelihood when the theory covariance is of the form given in Equation IA2I 
is: 

= Id^C-'CsC-'d - iTr (C-^Cb) (A3) 
where Tr is the trace operator. The curvature is: 

= -dtc-iCBC-^C^.C-M + \Tr (C-^CbC-^Cb,) (A4) 
oqBOq'B 2 

. It is straightforward to show that 

a^Ba = Tr (Baa^) (A5) 

for vector a and matrix B. One can use that identity to rewrite the first term in the curvature 
as follows: 

Tr [C-^CBC-^CB'C-^dd^] . (A6) 
If the correlation matrix is a good description of the data, we have {^ddF^ = C, or equivalently 



that C ^ddF ~ I/, the identity matrix. At this point, the standard treatment ( iBond et al. 
2 pool , e.g., ) is to replace C~^ddF with I which leaves 



^M£) ^ (C-'CbC-Cb.) . (A7) 

The most efficient implementation when using this approximation requires one to pre- 
calculate the set matrices C~^Cb, which requires an expensive matrix-matrix multiplication 
for each bandpower. One can then use the fact that Tr (AB) is an operation (essentially 
since one only need calculate the diagonal elements of AB). 

Instead of using l^ddF^ ^ C to cancel the first term in the curvature, MPILIKELY uses 
it to cancel the second term, leaving: 

The great advantage of this form of the approximate curvature is that it can be calculated us- 
ing strictly matrix-vector operations. In practice, on MPI clusters, we find that MPILIKELY 
typically spends 80-90% of its time inverting the covariance matrix, even when fitting dozens 
of parameters. In addition, storage requirements are halved since we don't have to store the 
set of matrices C^^C^. We also note that the error in Equation IA7I is exactly twice that 
of IA8I However, since the curvature only affects the path taken to the maximum likelihood 
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solution, both approximations will result in the same bandpowers. For final data products, 
we typically evaluate the full curvature at convergence for more accurate error bars. 

When finding the maximum-likelihood spect rum, MPILIKELY uses a modified version 



of the Levenberg-Marquardt algorithm (see e.g., iPress et al.l (119921 )) that generally reduces 
to Newton-Raphson iteration. The Levenberg-Marquardt control parameter A is initially set 
to zero, and remains there as long as the covariance remains positive-definite. On the first 
failure, it is set to unity, under the theory that if the current guess has overshot the maximum 
enough to make C non-positive-definite, a correction of the step by of order at least a factor 
of two is warranted. On continued failures, we increase A by a factor of 2, and on successes, 
we decrease it by \/2 until A falls below one, at which point we set it to zero. Our convergence 
criteria are that the largest step in any dimension, in terms of the error in that dimension, 
is less than some small fraction ( typically 0.01), and that A = during that iteration. In 
practice, A stays at zero, unless there is a power spectrum bin that is significantly negative. 
This happens when a random realization of the noise has substantially less power in it than 
expected, which pushes the power spectrum negative and introduces substantial skewness to 
the likelihood. 



A. 2. Spectrum From Matrices 

It is useful to be able to calculate the power spectrum and errors expected from either 
a noiseless data vector or a signal matrix. The treatment for the two cases is similar. We 
can invoke the identity in Equation IA51 replace ddF with the matrix D, and take advantage 
of the fact that Tr (AB) = Tr (BA) to rewrite the gradient as follows: 

= ^Tr (C-^DC-^Cb) - ^Tr (C-^Cb) . (A9) 

We can now find the expected spectrum from data drawn from covariance matrix D, marginal- 
ized over realizations of the data. When fitting to noiseless data, we have D = ddF + to 
marginalize over realizations of the noise. The gradient calculation is slowed by a factor of 
~ 3 since we now have to calculate C^^DC^^ instead of just C~^. The curvature calculation 
is slightly more complicated. Without an actual data vector, we cannot take advantage of the 
fast curvature calculation in Equation IA8I While one could use the curvature in Equation 
I A 71 this can become prohibitively slow in practice. MPILIKELY's solution is to draw a set 
of sample data vectors from D, and then average Equation IA8I for each of those realizations. 
The only operation this requires is a single initial Cholesky factorization of D to calculate 
the data vectors. 
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A.3. Source Projection 



The standard way of re moving; point sourc e s with known posit ions is by projecting them 



from the covariance matrix (IBond et al.l l2000l : iMason et al.l l2003l . e.g., ). This is roughly 



equivalent to masking out the source location in the map. In practice, this is done by adding 



Qsrc _ ^g^gT ^YiQ covariance matrix, where Si is the signal expected from the i source, 
and (3 is an (extremely large) amplitude. For sufficiently large /3, the power spectrum is 
insensitive to any signal proportional to Sj - i.e., we don't need to know the amplitude of 
the source. We have previously used large values of /3, the projection amplitude, but this 
can lead to numerical stability problems as large /5's cause C to become ill-conditioned. 
Instead, in MPILI KELY we take the analytic limit as /5 — oo using the Woodbury formula 



( jPress et al.lll992l . see e.g., ). For symmetric matrices, we have 



lim (C + /?SST) ^ = C-^-C~^S(S^C-^S) ^S^C^^ 



(AlO) 



One must also guard against the possibility of S being degenerate. If it is, the calculation 
of (S^C~^S) will fail. To deal with this possibility, we first scan S for repeated columns 
(same source entered twice). We then orthogonalize S through use of a QR factorization, 
and use the orthogonal matrix Q in Equation lAlOl instead of S. One could also use SVD, 
but QR is typically a factor of ~ 6 faster. With the addition of the QR factorization, the 
Woodbury formula never requires the inversion of an ill-conditioned matrix, and so as long 
as the source-free version of the covariance C is well-conditioned, source projection does not 
lead to numerical instability in calculating the source-projected C^^. 



A. 4. Likelihood Evaluations 

It is often useful to directly evaluate the likelihood of a given covariance matrix. Some 
instances where this is useful are when comparing the goodness-of-fit of two different models 
for the data, or when measuring quantities like confidence intervals that depend on the 
non-Gaussian nature of the likelihood surface. The analytic projection of sources using the 
Woodbury formula complicates likelihood evaluations because becomes singular (the 
projected modes are truly gone). The solution is to rotate into the subspace spanned by 
the complement of the source vectors. The rotation/compression matrix can also be found 
efficiently through a QR factorization of the source vectors (most QR implementations also 
have a way to find an orthogonalized matrix spanning the complement of the original matrix) . 
Each matrix (noise, banded signal...) and the data are then compressed. One could fit the 
spectrum with these compressed matrices, which gives the same spectrum one gets when 
using Equation lAlOl Since the matrices are smaller, the actual fitting of the power spectrum 
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is sped up. However, the initial overhead required (2 matrix-matrix multiphcations for each 
CMB band, plus typically a few others) is generally much larger than the total time to fit the 
spectrum with the uncompressed matrices, and so we in general only compress the matrices 
if we desire likelihoods. 



A. 5. Template Fitting 

MPILIKELY also supports simultaneous fitting of additive templates to the data when 

measuring the power spectrum. Wc use this for measuring the amplitude of foreground maps 
in the CBI data. For additive templates, the template model is subtracted from the data 
before calculating the likelihood, which then becomes: 

log(£) = -Vc-^d*-llog(|C|) (All) 

where d* — d — '^ajm^ for expected signal with amphtude aj and observed data d. 
Here, the index j runs over templates, not over data elements. The gradient and curvature 
of the CMB terms are unchanged, as long as they are calculated using d* instead of d. The 
template terms in the gradient are: 



^^°S^^)=m^tc-V. (A12) 



The curvature terms are: 



log (£) 
dajdaf 



(A13) 



and 



= -m^'C-CBC-'d. (A14) 

These curvature and gradient terms only require matrix-vector operations to calculate, and 
so have a negligible impact on the time required to fit the power spectrum. 



